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We report an analysis of charmless hadronic decays of charged B mesons to the final state K'^-k^ti^ , 
using a data sample of 470.9 ± 2.8 million BB events collected with the BABAR detector at the T{AS) 
resonance. We observe an excess of signal events, with a significance above 10 standard deviations 
including systematic uncertainties, and measure the branching fraction and CP asymmetry to be 
B {B+ K+n°TT°) = (16.2 ± 1.2 ± 1.5) x 10"^ and Acp {B+ K+tv^tt'-'') = -0.06 ± 0.06 ± 0.04, 
where the uncertainties are statistical and systematic, respectively. Additionally, we study the 
contributions of the B+ K* {S92) + tv° , B+ fo{9SO)K+, and B+ XcoK+ quasi-two-body de- 
cays. We report the world's best measurements of the branching fraction and CP asymmetry of the 
B+ K+Tv°TT° and B+ J^*(892) + 7r° channels. 

PACS numbers: 13.25.Hw, ll.30.Er, 12.39.-x 



I. INTRODUCTION 



Recent measurements of rates and asymmetries in 
B — >■ Kit decays have generated considerable interest be- 
cause of possible hints of new physics contributions [l|, Q ■ 
Unfortunately, hadronic uncertainties prevent a clear in- 
terpretation of these results in terms of physics beyond 
the Standard Model (SM). A data-driven approach in- 
volving measurements of all observables in the B — > Kn 
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system can in principle resolve the theoretical situation, 
but much more precise measurements are needed 0-11}. 

The ratios of tree-to-penguin amplitudes in the related 
pseudoscalar-vector decays B K*tt and B — > Kp are 
predicted to be two to three times larger than those 
in i? — > K-K. Hence, these decays could have consid- 
erably larger CP asymmetries and thus provide useful 
additional information jl-Q . In Table U we review the 
existing experimental measurements of the channels in 
the B — > K*TT system. Improved measurements of the 
K*+'!T^ decay can be obtained using the full r{AS) 
BABm. dataset. 



TABLE I: Experimental measurements of B — > K*Ti decays. 
Average values come from HFAG 10]. 



Mode 


B X 10" 


Acp 


References 


K*+%- 
K'+^° 
K'-°-K+ 
K*°-K° 


10.3 ± 1.1 
6.9 ±2.3 

^■^ -0.9 

2.4 ±0.7 


-0.23 ±0.08 
0.04 ± 0.29 ± 0.05 
-0 020 +° '"^'^ 

U.UZU _o.o61 

-0.15 ± 0.12 ± 0.02 


[15] 

[16, 17] 
[12,18] 
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The four K*'it decays populate six K-k-k Dahtz plots 
(the four Kp decays also produce four of the same six 
final states). To date, Dalitz plot analyses have been per- 
formed in the channels i^T+Tr+Tr" K2.tt'^it~ 
[ill and K+Ti--K° [H The first two of these have 
shown the presence of a poorly-understood structure, 
dubbed the /x(1300), in the tt+tt^ invariant mass dis- 
tribution. A study of the invariant mass spectrum in 
— > K^TT^TT^ decays could help elucidate the nature 
of this peak, since even-spin states will populate both 
/^vr+TT^ and Xtt'^tt*' (assuming isospin symmetry), while 
odd-spin states cannot decay to tt^tt'^. 



Knowledge of the dominant contributions to the 
i^T+TT^Tr" Dalitz plot may also help to clarify the inter- 
pretation of the inclusive time-dependent analyses [lil 
of — >■ K'^tt^tt^ [l^l- For such s penguin-dominated 
decays the naive Standard Model expectation is that the 
time-dependent CP violation parameter should be given 
by Sep ~ — ??cp sin(2/3), where ricp is the CP eigenvalue 
of the final state (+1 for -fC^Tr^Tr") and /3 is an angle of the 
Cabibbo-Kobayashi-Maskawa [Ij, [131 unitarity triangle. 
Currently, the results for — >■ i^gTr^Tr" show the largest 
deviation, among hadronic b ^ s penguin-dominated de- 
cays fld\, from the angle (3 measured in charmed decays, 
albeit with a large uncertainty. Such deviations could 
be caused by new physics, but in order to rule out the 
possibility of sizeable corrections to the Standard Model 
prediction, better understanding of the population of the 
Ktt'^tt'^ Dalitz plots is necessary. 



In this article, we present the results of a search for 
the three-body decay B^ K^tt'^tt^, including short- 
lived intermediate two-body modes that decay to this 
final state. A full amplitude analysis of the three-body 
decay would require detailed understanding of effects re- 
lated to the misreconstruction of signal events, such as 
the smearing of their Dalitz plot positions. These effects 
are significant in the final state under study, which in- 
volves two neutral pions. Therefore, in order to avoid 
heavy reliance on Monte Carlo (MC) simulations, we do 
not perform a Dalitz plot analysis, but instead extract 
information on intermediate modes including narrow res- 
onances (ii:*+(892)7r", /o(980)7^+ and XcoK+) by study- 
ing the two-body invariant mass distributions. 



There is no existing previous measurement of the 
three-body branching fraction, but several quasi-two- 
body modes that can decay to this final state have 
been seen, with varying significances. These include 
B+ f()i980)K+, observed in the /o(980) -> n+jT- 



channel [16|, [17| and also seen in /o(980) K^K 
B+ /2(1270)i^+, seen in /2(1270) ^ vr+Tr" [H, 
and B+ ^ K*+{892)tt°, seen in K*+{892) K+tt° 
The decay B~^ — » ~Xd)K^ has also been observed with 
and Xco ^ K+ R- 



XcO 



II. EVENT RECONSTRUCTION AND 
SELECTION 

The data used in the analysis were collected with the 
BaBAR detector [l^ at the PEP-II asymmetric-energy 
e+e^ collider at the SLAC National Accelerator Labo- 
ratory. The sample consists of an integrated luminos- 
ity of 429 fb~^ recorded at the T{AS) resonance ("on- 
peak") and 45fb~^ collected 40 MeV below the resonance 
("off-peak"). The on-peak data sample contains the full 
BABm T(45) dataset, consisting of 470.9 ± 2.8 million 
BB events. 

We reconstruct B'^ — >■ K'^tt'^tt'^ decay candidates by 
combining a candidate with two neutral pion can- 
didates. The K'^ candidate is a charged track with 
transverse momentum above 0.05 GeV/c that is consis- 
tent with having originated at the interaction region. 
Separation of charged kaons from charged pions is accom- 
plished with energy-loss information from the tracking 
subdetectors and with the Cherenkov angle and num- 
ber of photons measured by a ring-imaging Cherenkov 
detector. The efficiency for kaon selection is approxi- 
mately 80 % including geometrical acceptance, while the 
probability of misidentification of pions as kaons is below 
5 % up to a laboratory momentum of 4 GeV/ c. Neutral 
pion candidates are formed from pairs of neutral clus- 
ters with laboratory energies above 0.05 GeV and lat- 
eral moments [1^ between 0.01 and 0.6. We require 
the mass of the reconstructed 7r° to be within the range 
0.115 GeV/ < nij^ < 0.150 GeV/ and the absolute 
value of the cosine of the decay angle in the 7r° rest frame 
to be less than 0.9. Figure [T] shows the distribution of 
the mass of neutral pion candidates in on-peak data. Fol- 
lowing this selection, when forming the B candidate, the 
tt" candidates have their masses constrained to the world 
average value ^2^. 



^ 40000 
> 

S 30000 



20000 



g 10000 
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0.12 0.14 0.16 
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FIG. 1: Masses of tv^ candidates in on-peak data. The arrows 
indicate the selection requirements. 

We exclude candidates consistent with the B~^ ^• 
K'^K^, K'^ tt'^tt'^ decay chain by rejecting events with 
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a pair of tt" mesons that satisfies 0.40 GeV/c^ < 771^07^0 < 
0.55 GeV/c^. This veto has a signal efficiency of at least 
96 % for any charmless resonant decay and of almost 
100 % for nonresonant B+ K+n°n° and B+ XcoK+ 
decays. 

Because of the presence of two neutral pions in the final 
state, there is a significant probability for signal events 
to be misreconstructed, due to low momentum photons 
that are replaced by photons from the decay of the other 
B meson in the event. We refer to these as "self-cross- 
feed" (SCF) events, as opposed to correctly reconstructed 
(CR) events. Using a classification based on Monte Carlo 
information, we find that in simulated events the SCF 
fraction depends strongly on the resonant substructure 
of the signal, and ranges from 2 % for 5+ — ^ XcoK^ 
decays to 30% for B+ /2(1270)if+ decays. 

In order to suppress the contribution arising from 
the dominant background, due to continuum e^e~ — > 
qq [q — u, d, s, c) events, we employ a neural network 
that combines four variables commonly used to discrimi- 
nate jet-like qq events from the more spherical BB events. 
The first of these is the ratio of the second-to-zeroth order 
momentum- weighted Legendre polynomial moments, 

E H3cos2 0,-l)p, 

^2 _ igROE /,N 

ieROE 

where the summations are over all tracks and neutral 
clusters in the event excluding those that form the B 
candidate (the "rest of the event" or ROE) , pi is the par- 
ticle momentum and 9i is the angle between the particle 
and the thrust axis of the B candidate. The three other 
variables entering the neural network are the absolute 
value of the cosine of the angle between the B direction 
and the beam axis, the absolute value of the cosine of 
the angle between the B thrust axis and the beam axis, 
and the absolute value of the output of a neural network 
used for "flavor tagging" , i.e. , for distinguishing B from 
B decays using inclusive properties of the decay of the 
other B meson in the T{4S) BB event [sO]. The 
first three quantities are calculated in the center-of-mass 
(CM) frame. The neural network is trained on a sam- 
ple of signal MC and off-peak data. We apply a loose 
criterion on the neural network output (NNout), which 
retains approximately 90 % of the signal while rejecting 
approximately 82 % of the qq background. 

In addition to NNout , we distinguish signal from back- 
ground events using two kinematic variables: 

mES = ^Ei-pl, (2) 

AE = E*B-V^/2, (3) 

where 

Ex = {s/2 + p,-pb)/E,, (4) 

y/s is the total CM energy, {Ei, p^) and (i^s, Ps) are the 
four-momenta of the initial e'^e~ system and B candi- 
date, respectively, both measured in the lab frame, while 



the star indicates the CM frame. The signal toes dis- 
tribution for CR events is approximately independent of 
the i?+ — > +7r°7r° Dalitz plot distribution and peaks 
near the B mass with a resolution of about 3MeV/c^. 
We select signal candidates with 5.260 GeV/c^ < toes < 
5.286 GeV/c^. The CR signal AE distribution peaks 
near zero, but has a resolution that depends on the 
event-by-event Dalitz plot position, the PDF of which 
is a priori unknown. Prior to the selection of multi- 
ple candidates (see below), we make the requirement 
\AE\ < 0.30 GeV, in order to retain sidebands for back- 
ground studies. However, to avoid possible biases [Slj 
we do not use AE' in the fit described below and instead 
apply tighter selection criteria for events entering the fit, 
-0.15 GeV < AE < 0.05 GeV. These criteria have an 
efficiency of about 80 % for signal while retaining only 
about 30 % of the background, both compared to the 
looser requirement |Ai?| < 0.30 GeV. 

The efficiency for signal events to pass all the selection 
criteria is determined as a function of position in the 
Dalitz plot. Using an MC simulation in which events 
uniformly populate phase space, we obtain an average 
efficiency of approximately 16 %, though values as low as 
8% are found near the corners of the Dalitz plot, where 
one of the particles is soft. 

An average of 1.3 B candidates is found per selected 
event. In events with multiple candidates we choose the 
one with the smallest value of a variable formed from 
the sum of the values of the two tt*^ candidate masses, 
calculated from the difference between the reconstructed 
TT^ mass with respect to the nominal tt^ mass. This pro- 
cedure has been found to select the best reconstructed 
candidate more than 90 % of the time, and does not bias 
our fit variables. 

We study residual background contributions from BB 
events using MC simulations. We divide these events 
into four categories based on their shapes in the TOes 
and AE distributions. The first category comprises two- 
body modes (mainly i?+ iiT+Tr"); the second contains 
three-body modes (mainly B~^ — )■ K*'^{—>- K~^Tr'^)'j and 
5+ — >■ TT+TT^Tr"); the third and fourth are composed of 
higher multiplicity decays (many possible sources with or 
without intermediate charmed states) with missing par- 
ticles, and are distinguished by the absence or presence 
of a peak in the TOes distribution, respectively. Based on 
the MC-derived efficiencies, total number of BB events, 
and known branching fractions [l^jll^j we expect 70 ±9, 
39 ± 18, 1090 ± 40 and 170 ± 30 events in the four cate- 
gories respectively. 



III. STUDY OF THE INCLUSIVE B+ K+TT°ir° 
DECAY 

To obtain the B'^ — >■ K'^tt'^tt'^ signal yield, we perform 
an unbinned extended maximum likelihood fit to the can- 
didate events using two input variables: toes and NNout- 
For each component j (signal, qq background, and the 
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four BB background categories) , we define a probability 
density function (PDF) 



P] = 7',(mEs*)7',(NNout'), 



(5) 



where the index i runs over the selected events. The 
signal component is further separated into CR and SCF 
parts 



= (1 - /sCF)^CR(™Es')^CR(NNo,t') + 

/scfT^scf (mEs* )7^scf (NN 

out ) 7 



(6) 



where /scf is the SCF fraction. The extended likelihood 
function is 



(7) 



where rij^k) is the yield of the event category j (fc). 

For the signal, the tties PDFs for CR and SCF are de- 
scribed by an asymmetric Gaussian with power-law tails 
and a third-order Chebyshev polynomial, respectively. 
Both CR and SCF NNout PDFs are described by non- 
parametric PDFs (one-dimensional histograms). We fix 
the shape parameters of the signal rriBs PDFs to the val- 
ues obtained from the i?+ — > i^+Tr^Tr" phase-space MC 
sample. The parameters are corrected to account for pos- 
sible differences between data and MC simulations, us- 
ing correction factors determined with a high-statistics 
control sample of B+ D^p+ {K+n-TT°) (tt+tt") 
decays. For the continuum background, we use an AR- 
GUS function [s^ to parameterize the tties shape. The 
endpoint of the ARGUS function is fixed to 5.289 GeV/c^ 
whereas the shape parameter is allowed to float in the 
fit. The continuum NNout shape is modelled with a 
20 bin parametric step function, i.e. , a histogram with 
non- uniform bin width and variable bin content. One- 
dimensional histograms are used as nonparametric PDFs 
to represent all fit variables for the four BB background 
components. The free parameters of our fit are the yields 
of signal and continuum background together with the 
parameters of the continuum toes and NNout PDFs. All 
yields and PDF shapes of the four BB background cat- 
egories are fixed to values based on MC simulations. 

The results of the fit are highly sensitive to the value of 
/scf, which depends strongly on the Dalitz plot distribu- 
tion of signal events and cannot be determined directly 
from the fit. To circumvent this problem, we adopt an 
iterative procedure. We perform a fit with /scf fixed 
to an initial value. We then construct the signal Dalitz 
plot from the signal probabilities for each candidate event 
[sWeights) calculated with the /Plot technique [S^l, and 
determine the corresponding average value of /scf. We 
then fit again with /scf fixed to the new value, and re- 
peat until the obtained values of the total signal yield 
(CR -I- SCF) and /scf are unchanged between iterations. 
This method was validated using MC and was found to 
return values of /scf that are accurate to within 3 % of 



the nominal SCF fraction. Convergence is typically ob- 
tained within three iterations. 

We cross-check our analysis procedure using the high 
statistics control sample described above. We impose se- 
lection requirements on the D and p candidates' invari- 
ant masses: 1.84GeV/c^ < m^+^-jrO < 1.88GeV/c^ and 
0.65GeV/c2 < m^+^o < 0.85 GeV/c^ We fit the on-peak 
data with a likelihood function that includes components 
for the control sample, all BB backgrounds, and qq. We 
find a yield that is consistent with expectation based on 
the world-average branching fractions 29] . 

We apply the fit method described above to the 31 673 
selected candidate B~^ — > K'^tt'^'tt^ events. Convergence 
is obtained after four iterations with a yield of 1220 ± 85 
signal events and a SCF fraction of 9.7%. The results of 
the fit are shown in Fig. O The statistical significance of 
the signal yield, given by \/2Aln£ where Aln£ is the 
difference between the negative log likelihood obtained 
assuming zero signal events and that at its minimum, is 
15.6 standard deviations (cr). Including systematic un- 
certainties (discussed below), the significance is above 
10 cr. 

To obtain the — > A'+7r°7r° branching fraction us- 
ing the result of the fit, we form, for each event, the ra- 
tio of the signal ^Weig/it and the efficiency determined 
from its Dalitz-plot position. Summing these ratios over 
all events in the data sample, we obtain an efficiency- 
corrected signal yield of 7427 ± 518 events. The .^yVeiqht 
calculation accounts for the fixed BB backgrounds [33|. 
The Dalitz plot distributions obtained before and after 
applying the efficiency correction are shown in Fig. |31 
We apply further corrections for the effect of the K'^ veto 
(98%); differences between data and MC for the 7r° re- 
construction efficiency, determined from control samples 
of T decays as a function of tt'^ momentum (95.7%); and 
a bias in the fitted signal yield (raw bias 44 events), as 
determined from Monte Carlo pseudoexperiments gener- 
ated with a signal component with the same values of the 
yield and SCF fraction as found in the fit to data. Fi- 
nally, we divide by the total number of BB events in the 
data sample to obtain our measurement of the branching 
fraction B {B+ K+tt^tt^) = (16.2 ± 1.2 ± 1.5) x 10-^ 
where the first uncertainty is statistical and the second 
is systematic. 

The systematic uncertainty includes contributions 
from the PDF shapes; the fixed BB background yields; 
the estimation of the SCF fraction; intrinsic fit bias; se- 
lection requirements; and the number of BB pairs in 
the data sample. Here we provide further details on 
each of these sources of systematic uncertainty and de- 
scribe briefly how each is evaluated. A combined uncer- 
tainty for CR signal and BB background NNout PDF 
shapes (4.9%) is evaluated using uncertainties in the 
data/MC ratio determined from the S"*" D^p^ — 
(if+TT^Tr") (tt+tt*^) control sample and applying them 
simultaneously to the CR signal and BB background 
NNout PDFs. The same control sample is used to eval- 
uate the uncertainties in CR signal mss PDF shapes 



5:26 5.27 5.28 

mj,s (GeV/c^) 

FIG. 2: Projections of candidate events onto tues (left) and NNout (right), following requirements on the other fit variable that 
enhance signal visibility. These requirements retain 60% of signal events for the tues plot and 87% of events for the NNout plot. 
Points with error bars show the data, the solid (blue) lines the total fit result, the dashed (green) lines the total background 
contribution, and the dotted (red) lines the qq component. The dash-dotted lines represent the signal contribution. 
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FIG. 3: Signal Dalitz plot distributions obtained using sWeights before (left) and after (right) efficiency correction. In order to 
define a unique position for each event, the Dalitz plot is shown as '^^^+^o vs. m^o^o , where JTI^+^o is the smaller of the two 

if'^Tr" invariant masses. Resonance bands are visible for AT* (892)"*" at m^jf+^o ~ 0.8 GeV^/c*, /o(980) at m^o^o ~ 1 GeV^/c'* 

and Xco at m^o^o ~ 12 GeV^/c"*. 



(0.8%). The uncertainty in the SCF fraction (2.5%) is 
estimated by varying the value used in the fit within a 
range of uncertainty determined from Monte Carlo pseu- 
doexperiment tests of our iterative fitting procedure. Un- 
certainties in the SCF signal itibs and NNout PDF shapes 
(1.7% and 0.7%, respectively) are evaluated by consid- 
ering a range of SCF shapes corresponding to different 
signal Dalitz plot distributions. An uncertainty in the 
correction due to fit bias (1.9%) is assigned, which cor- 
responds to half the correction combined in quadrature 
with its error. Uncertainties in the BB background toes 
PDF shapes due to data/MC differences (1.6%) are eval- 
uated by smearing the PDFs with a Gaussian with pa- 
rameters determined from the B~^ — >■ control sam- 
ple. The uncertainties in the BB background PDFs due 



to finite MC statistics (0.8 %) are determined by vary- 
ing the contents of the bins of the histograms used to 
describe the PDFs within their errors. Uncertainties in 
the fixed BB background yields (1.4%) are evaluated 
by varying these yields within their uncertainties. Con- 
tributions to the uncertainty in the selection efficiency 
arise from the /\E (4.0 %) and NNout (3.0 %) selection re- 
quirements, neutral pion reconstruction (2.8%), the 
veto correction (2.0%), kaon identification (1.0%) and 
tracking (0.4%); The uncertainty in the number of BB 
pairs in the data sample is 0.6 %. Including only system- 
atic uncertainties that affect the fitted yield, the total is 
6.5%. The total systematic uncertainty on the branch- 
ing fraction is 9.0 %. Table HIl summarizes the systematic 
contributions. 
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TABLE II: Summary of systematic uncertainties for the in- 
clusive branching fraction measurement. 



Source 


Uncertainty 


CR signal and BB background NNout PDFs 


4.9% 


CR signal tties PDF 


0.8% 


SCF fraction 


2.5 % 


SCF signal tties PDF 


1.7% 


SCF signal NNom PDF 


0.7% 


Fit bias 


f.9% 


BB background rriEs PDFs 


1.6% 


BB background PDFs (MC statistics) 


0.8% 


BB background yields 


1.4% 


Subtotal 


6.5% 


AE selection efficiency 


4.0% 


NNout selection efficiency 


3.0% 


Neutral pion efficiency 


2.8% 


K'^ veto 


2.0% 


Particle identification efficiency 


1.0% 


Tracking efficiency 


0.4% 


Nbb 


0.6% 


Total 


9.0% 



The CP asymmetry is measured as 

Nb- - Nb+ 



A, 



CP 



(8) 



where Nb+(b-) is the nmnber of events from 
5+ — )■ if +7r°7r° (CP conjugate decay) and is obtained 
by including in the above-described fit the value of the 
kaon charge. The fit returns an asymmetry of Acp — 
—0.06 ±0.06 ±0.04. Most of the systematic uncertainties 
that affect the branching fraction cancel in the asymme- 
try. However, the following sources are considered and 
evaluated for the Acp measurement. Detector-induced 
asymmetries have been studied in previous similar anal- 
yses [IB 13 ^^'^ found to be small (0.5%). We evaluate 
the possibility that our selection induces an asymmetry 
by measuring the CP asymmetry in the B'^ — )• D'^p^ 
control sample (3.0 %), where none is expected. The BB 
background asymmetries are fixed in our fit; the uncer- 
tainty from this is evaluated (1.8%) by varying these by 
a weighted average of the CP asymmetries of the con- 
tributing BB decays. Finally the fit bias is estimated 
from MC pseudoexperiments (1.2%). 



around 1550MeV/c^ that could arise from interference 
between various resonances in this region. A broad peak 
around 1400MeV/c^ in the K^tt'^ invariant mass distri- 
bution could be due to contributions from spin-0 and/or 
spin-2 if* (1430)+ states. 



The numbers of signal events for the quasi-two-body 
contributions are determined by defining signal regions 
around the peaks of the resonances. Efficiency-corrected 
sWeights are summed in the same way as used to mea- 
sure the inclusive branching fraction. To estimate contri- 
butions from nonresonant and resonant B'^ — )■ if +7r''7r'' 
decays other than the quasi-two-body decays under con- 
sideration (which we refer to as background in this sec- 
tion), the same procedure is applied to sidebands on ei- 
ther side of each signal region in the two-particle invari- 
ant mass. The background yields are estimated as the 
normalized averages of the two sidebands' yields and are 
subtracted from the efficiency-corrected yields in the sig- 
nal regions. The signal and sideband regions are illus- 
trated by arrows for each of the three quasi-two-body 
modes in Fig. [5] We use this approach rather than a full 
Dalitz plot analysis since the latter would require more 
detailed understanding of the properties of SCF events. 
Our method does, however, suffer from systematic un- 
certainties (evaluated below) due to other contributions 
to the Dalitz plot and possible interference effects. This 
precludes its use for studying quasi-two-body decays via 
broad resonances. We have validated our approach using 
ensembles of MC simulations with varying mixtures of 
resonant substructure, and find that in all cases we are 
able to correctly obtain the true values of the branch- 
ing fractions of the quasi-two-body decays under study, 
which all have narrow intermediate states under study. 



Fits to the efficiency-corrected invariant mass distri- 
butions are used to cross-check the results of the sub- 
traction method. In these fits we describe the signal dis- 
tributions with double-Gaussian functions, with parame- 
ters obtained from MC simulations, and the background 
shapes with polynomials. The two methods yield consis- 
tent results, both in MC simulations and in data. 



IV. STUDY OF QUASI-TWO-BODY 
CONTRIBUTIONS 

We use the sPlot distributions obtained from the fit 
and projected onto the Dalitz plot axes to search for 
peaks from intermediate resonances. These projections 
are shown for both if +7r*' and tt^tt*^ invariant masses in 
Fig.H Signal peaks from if*(892)+, /o(980) and Xco are 
clearly observed. We do not see any enhancement that 
could be attributed to the /x(1300), though the tt^tt" 
invariant mass distribution contains a pronounced dip 



After background subtraction we obtain 
efficiency-corrected signal yields of 1078 ± 197 for 
B+ if*(892)+7r°, 1186 ± 241 for B+ /o(980)if+, 
and 245 ± 105 for i?+ — )• XcoK'^. We correct each 
yield for the inefficiency of the corresponding signal 
region selection, obtained from Monte Carlo simulations. 
Finally the yields are corrected: (i) for bias, estimated 
from Monte Carlo pseudoexperiments; (ii) for tt*^ effi- 
ciency, using the momentum distributions of both tt'^ 
mesons from a Monte Carlo cocktail refiecting the yields 
obtained in data; (iii) in the case of the if* (892)+ yield 
only, for the if^ veto. Finally, we divide by the number 
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FIG. 4: Signal sP^ot distributions not corrected for efficiency for (a) 0.5 GeV/c^ < m^+^o . < 2.0GeV/c^, (b) 0.5 GeV/c^ < 
m^o^o < 2.0GeV/c^ and (c) 3.0GeV/c^ < m^o^o < 4.0GeV/c^. m^+^o^.^ is the K^n'^ combination witii lower invariant mass. 
Excesses of events in the /o(980), XcO, if*(892)+, and A'*(1430)+ mass regions are clearly visible. 



of BB pairs to obtain the product branching fractions 

B {B+ i^*(892)+7rO) x B (X*(892)+ ^ K+tt^) = 

(2.7 ±0.5 ±0.4) X 10-^ 
B {B+ /o(980)if +) X B (/o(980) ^ tt^tt") = 

(2.8 ±0.6 ±0.5) X 10-6, 
B {B+ ^ Xc^K+) X B (xco -> ^°^°) = 

(0.51 ±0.22 ±0.09) X 10-6, 

(9) 

where the first uncertainties are statistical and the sec- 
ond systematic. The sum of these contributions does not 
saturate the inchisive branching fraction, indicating sig- 
nificant contributions from other sources, as is also clear 
from Fig. [3] and Fig. IH and expected from the results of 
studies of B+ -> K+tt+tt" decays [l6l[l7j. 

Systematic uncertainties include all the same sources 
in the same relative amounts as evaluated for the inclu- 
sive decay except for fit bias, veto and 7r° efficiency, 
which are evaluated separately for each quasi- two- body 
mode. We also evaluate the following additional contri- 
butions. The uncertainty due to the method of back- 
ground subtraction (3.5% for i^*(892)+7rO, 11.9% for 
/o(980)X+ and 13.5% for XcnK+) is obtained by com- 
paring the nominal results with those from obtained with 
alternative sideband regions. We evaluate the potential 
effect of interference (10.0%, for /o(980)i4'"'' only) using 
toy Monte Carlo events generated for a Dalitz plot model 
containing /o(980) and nonresonant components with rel- 
ative magnitudes obtained from the fit results, and a rel- 
ative phase sampled in a range that gives distributions 
consistent with the data. Finally we consider possible 
data/MC differences affecting the signal region efficiency 
correction (5.6% for is:*(892)+7r°, 3.8% for /o(980)ii'+, 
and 0.4% for XcoK^) determined from the change in the 
result when the SCF fraction is varied in Monte Carlo 
events. The X*(892)+7r° and XcoK^ branching fraction 
measurements are not affected by systematics due to in- 
terference. For the former, effects of interference with 
K^TT^ S-wave contributions cancel when integrated over 
the part of the Dalitz plot inside the signal mass win- 
dow, while P-wave contribution are not expected based 
on studies of related decays [l^ [l3|. For the latter, the 
small width implies that interference will be negligible. 



TABLE III: Summary of systematic uncertainties for the 
branching fraction measurement of the quasi-two-body res- 
onances. The breakdown of the systematics affecting the in- 
clusive branching fraction measurement is given in Table HIl 



Source 


Uncertainty (%) 






if* (892)-' 


-n° fo{9S0)K+ 




Subtotal from inclusive 


8.1 


8.1 


8.1 


Background subtraction 


3.5 


11.9 


13.5 


Interference 




10.0 




Fit bias 


6.6 


2.1 


6.8 


Mass cut efficiency 


5.6 


3.8 


0.4 


Tv'^ efficiency 


3.1 


3.5 


2.6 


Ks veto 


2.0 






Total 


12.9 


18.4 


17.4 



A list of the systematic uncertainty contributions is given 
in Table Uni 

To obtain the B decay branching fractions, we correct 
for B {K*{892)+ K+tt°) = 1/3 and B {xcO t^°t^°) = 
(8.4 ± 0.4) X 10-3 X 1/3 [291], where the factors of 1/3 are 
due to isospin. (The branching fraction of /o(980) — > 
TT^TT*^ is unknown, hence we cannot correct for it.) We 
obtain 

B{B+ ^ K* (892)+ n^) = (8.2 ± 1.5 ± 1.1) x lO'^ , 
B{B+^XcoK+) = (18±8±3± 1) X 10-^ 

(10) 

where the first uncertainty is statistical, the second sys- 
tematic, and the third (for 5+ XcoK~^) is from the 
subdecay branching fraction. 

We obtain the CP asymmetries of the quasi-two-body 
modes with the same method used to obtain the quasi- 
two-body branching fractions, except we distinguish the 
yields of the i?+ and i?- decays. We obtain the following 
asymmetries: 

Acp {B+ K*{892)+n") = -0.06 ± 0.24 ± 0.04 , 
Acp {B+ -> /o(980)/-i:+) = 0.18 ± 0.18 ± 0.04 , 
Acp {B+ XcoK+) = -0.96 ± 0.37 ± 0.04 , 

(11) 

where the first uncertainty is statistical and the sec- 
ond systematic. The sources of systematic uncertainty 
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FIG. 5: Efficiency-corrected signal (vertical red arrows) and sideband (horizontal blue arrows) regions around (a) K* (892)'^ , 
(b) /o(980) and (c) XcO invariant mass peaks. The curves show the results of the fit used to cross-check the procedure, for the 
total (blue continuous line) and background-only (dashed red line) components. 



are the same as for the inclusive CP asymmetry mea- 
surement. The measurements of CP asymmetries for 
foi980)K+ and B+ XcoK+ are consistent with 
the world average values based on decays of the inter- 
mediate resonances to tt+tt" [i3, 29]. The B'^ — >■ XcoK^ 
result has a large and non-Gaussian uncertainty and its 
difference from zero is not statistically significant. 

V. CONCLUSION 

In summary, using the full BaBAR data sample of 
429 fb"^ collected at the T(AS) resonance, we observe 
charmless hadronic decays of charged B mesons to the 
final state K'^tt^tt^ . The signal has a significance above 
10 cr after taking systematic effects into account. 

We study the Dalitz plot distribution of the signal 
events, and do not see any excess that could be attributed 
to the /x(1300). However, due to the possible compli- 
cated interference pattern, we cannot draw any strong 
conclusion about this state from our analysis. We mea- 
sure the product branching fractions and direct CP asym- 
metry parameters of the quasi-two-body modes with nar- 
row resonance peaks in the K'^tt^tt^ Dalitz plot. 

The results are summarized in Table ITVl All measured 
CP asymmetries are consistent with zero. The branch- 
ing fraction result for — >• XcoK^ is consistent with 



the world average, while that for B+ if* (892) +71° is 
consistent with and more precise than our previous mea- 
surement [iBl, which it supersedes. 
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